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Abstract 



O 
O 

£N) ' We have studied in detail the spectrum of the decay \i — > e v e v^ as a function of the electron 

helicity, and verify the prediction Q that a significant fraction of the electrons in this decay is right- 
j^' 1 handed. These "wrong-helicity" electrons persist in the limit A = m e /m M — » 0, and are connected 

with helicity- flip bremsstrahlung in QED. The longitudinal polarization of the electron is calculated as 
a function of the photon and electron energy, and deviates systematically from the naive V-A prediction 
Pl = — 1. The right-handed component is concentrated in the collinear region 6 < m e /E e . In the limit 
A — > 0, we reproduce the results obtained in Q using the helicity-flip splitting function introduced by 
Falk and Sehgal. 

> . 

m . 1 Introduction 
(N 
O 

' radiative decay [i~ — > e~D e u tt 'y are not purely left-handed, even in the limit m e — > 0. The decay width into 

■^j- , right-handed electrons was shown to be Tr = j^^o, where Tq = G|m^/(1927r 3 ). The energy spectrum of 

these wrong-helicity electrons was calculated, as also that of the photons accompanying them. These results 
, were obtained using the helicity-flip splitting function Dhf(z) = f^z introduced by Falk and Sehgal 

The appearance of right-handed electrons in ^t-decay is, at first sight, surprising since it goes against 
the conventional wisdom based on the V-A structure of weak interaction and the common assumption of 
hclicity-conservation in massless QED. However persistence of helicity-flip bremsstrahlung in the m e — > 
limit has been established in several calculations of radiative processes |S1 E] El HO 0, going back to the 
seminal paper of Lee and Nauenberg jH] . In all cases, the origin of the effect is the characteristic behaviour 
of helicity-flip bremsstrahlung at small angles 



In a recent Letter £Q it was pointed out that as a consequence of helicity-flip bremsstrahlung, electrons in the 



/m e \ 2 d^ 

S3: d ^bd [^ + (if) 2 f 

leading to a non-zero integrated rate that survives the limit m e — > 0. 

In the present paper, we derive the hclicity-dcpcndcnt spectrum of the decay pT — > e~i> e v^ from first 
principles, without neglecting the electron mass. We obtain the longitudinal polarization Pl of the electron 
as a function of the electron energy, photon energy and the angle between the electron and photon. We then 
consider the limit m e — > to demonstrate that the deviation from Pl = — 1 is present even in the chiral 
limit of QED. The results are compared with those obtained in the equivalent particle approach based on 
the helicity-flip fragmentation function, thus providing a confirmation of the statements in Ref. 1 . 



2 Matrix Element and Phase Space for fi — > e v e v^ 

The matrix element for the decay ji~ — > e~P e i^7 may be obtained from the Feynman diagrams shown in 
Fig-GJ- Adopting the notation of Ref. [51 ITU] . we obtain the square of the matrix element, summed over the 
photon polarization, taking the muon to be unpolarized. The spectrum is derived in terms of the invariants 

2pf,-p e 2p^-k 2p e -k 

x := ^ — ' V := T~ ' z := T~ ( 2 ) 

mf. mf. mf. 
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Figure 1: Feynman diagrams 



and the mass ratio A := m e /m^. In the rest frame of the muon, 



a; = — , y = ^ , z = ^ b - cos(#) V^ 2 - 4A 2 ] , (3) 
m M m M 2 

i^e, E 1 being the energies of the electron, photon and cos(0) the angle between them. The boundary of phase 
space is determined by the condition —1 < cos(0) < 1, and < Q 2 < (1 — A) 2 m 2 , Q 2 being the "missing 
mass" of the neutrino pair. This leads to the condition 

o<y<ym(x), 2\< x <i + \ 2 1 (4) 

where 

Vm(x) = 2 1 + A " = X . (5) 

V ' 2-a; + Va; 2 -4A 2 V ; 

Equivalently, the limits of phase space may be expressed as 

2A< x<x m (y), 0<y<l-A, (6) 

where 

A 2 + (l-y) 2 

x m (y) = . (7) 

1-2/ 

The differential decay width for pT — ► e~P e u fl 'y in the variables x, y and cos((9) is then given by 

U , ^ ff) . ) = (-1) r Vx 2 - 4 A 2 — 2 \ (.go ± 9x) (8) 
Vda;dj/dcos(y)/R,L 2 7r yz z 2 

where R and L denote right-handed and left-handed electrons, corresponding to the polarization vectors 

, s /|Pe| E e \ /|p e | E e „ \ 

( s e) R = , Pe Se L = , Pel- (9) 

The functions go and <?i are given by 

go := — 2 x y z 2 + 6 x y z 3 — 6 x y 2 z — 8 x y 2 z 2 + 6 x y 3 z + 6 x z 2 + 8x z 3 — 6 x 2 y z ~ 8 x 2 y z 2 

+ 8 x 2 y 2 z - 4 x 2 z 2 + 4 a; 3 y z + 6 y z 2 + 5 y z 3 - 2 y z 4 - 2 y 2 z 2 + 2 y 2 z 3 - 3 y 3 z - 2 y 3 z 2 
+ 2 y 4 z - 6 z 3 - 4 z 4 

+ A 2 (8xyz + 6xyz 2 + 2iry 2 z-|-6;ry 2 — 8 a; y 3 + 6 x z 2 — 6 a; 2 y z — 4 a; 2 y 2 + 6 y z 2 — 3 y z 3 

- 6 y 2 z - 2 y 2 z 2 + 5 y 3 z + 6 y 3 - 4 y 4 - 8 z 2 - 6 z 3 ) 
+ A 4 (6.xy 2 -6y 2 z-8y 2 + 6y 3 ), 
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<7x := {6a; yz — 4 a; yz + 6x y z — 8xy z+3xyz—6x y z—2xy z — 6x z +Ax z 
-~6xyz 2 + 2x 2 yz 2 + 8x 3 yz 2 + 2xy 2 z 2 + 8x 2 y 2 z 2 + 2xy 3 z 2 + 6xz 3 -8x 2 z 3 

o 2 3 2 3 4 4 

— hxyz —6a; yz — 2xy z +Axz +2xyz 

+ X 2 (-6 x 2 y 2 + 4 x 3 y 2 - 6 x y 3 + 8 x 2 y 3 - 6 y 4 + 8 x y 4 - 24 x y z + 16 x 2 y z ~ 4 y 2 z 3 

+ 4 y 5 + 2 x 3 y z + 22 x y 2 z - 6 x 2 y 2 z + 4 y 3 z - 7 x y 3 z - 4 y 4 z + 24 z 2 - 16 x z 2 - 2 x 2 z 2 
- 16 y z 2 - 18 x y z 2 - 2 x 2 y z 2 - 6 y 2 z 2 + 10 x y 2 z 2 + 4 y 3 z 2 + 16 z 3 + 2 x z 3 + x y z 3 ) 
+ A 4 (24 y 2 - 16. xy 2 - 2x 2 y 2 - 16 y 3 - 2xy 3 - 2y 4 - 8xyz + 16y 2 z + 2xy 2 z - Ay 3 z 
+ 8z 2 -2y 2 z 2 ) 
1 



+ 8AV} 



2 \lx 2 - 4 A 2 ' 

(11) 



Adding the expressions for dr^ and dlY, we obtain the spectrum for unpolarized electrons, which depends 
on the function g only |1U|. The specific effects of electron helicity are contained in the function g\. 

3 Unpolarized Spectra (summed over electron helicity) 

Summing over the electron helicity, the differential decay rate of — > e~v e y i /y is 

[ A , ^ m ) =(-l)r ^V^ 2 ~Ty^g . (12) 
Vda;dydcos(#)/R+L 2ir y z 2 

From this we can derive the following spectra. 

3.1 Dalitz plot in electron and photon energies 

(-^r) =(-l)r ^- f +1 d cos(d) ^x 2 - 4 A^ -L go (x, y, z; A) 
Vdxdy/R+L 2 ir J_ 1 yz 2 

= r -^- (.4{l2y(-6 + 3y + y 2 )+:r(-72 + 78y + 33y 2 -6y 3 ) 
D7T y V 



y 

+ 2x 2 (24+12y-5y 2 + 2y 3 ) 



+ A 2 [96 - 72y + 4y 2 - 4 y 3 + 9 a; (-8 - 2 y + y 2 )]} 
+ 12 C { - 4 a; 3 + x 2 (6 - 8 y) + 6 A 4 y - 6 x (-1 + y) y + (3 - 2 y) y 2 
+ A 2 [6a; 2 + (6-5y)y-2a;(4 + y)]}) , 



(13) 



where we have defined 



A := Vx 2 -4X 2 , £:=-log 



2 °Va;-„4 



(14) 



3.2 Photon Energy Spectrum 

The photon energy spectrum (dr/dy)R+L is given explicitly in Ea. l)A.l|) of the Appendix The leading term 
for small y (small photon energies) is 

/dT\ J^r — - r(-17 + 64A 2 -64A 6 + 17A 8 ) + (-12 + 144A 4 -12A 8 ) log(A)l . (15) 
Vdy/R+L 6n y L J 
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In the limit of small A, this reduces to the approximate form 



/dr 



v . / — r JL 1 [ - 17 - 12 log(A)] 
V dy / r+l 6 7r y 



(16) 



which agrees with the old calculation of Eckstein and Pratt |TTJ and Kinoshita and Sirlin J5| > m which terms 
of order A 2 and higher powers were neglected. The integrated rate for y > y is 



, .> \->0,yo->0 „ a 

^R+h(yo; A) > r - 

7T 



loe 



g(A)(2 log(yo) + lo g'M - y 



13409 



2520 



(17) 



which agrees with Ref. except in the constant (non- logarithmic) term, which was — 7r 2 /6 + 601/144 in 
the old calculation. 



3.3 Electron energy spectrum 

The exact result for the electron energy spectrum (dr/dx)R+L is given in Ea. (|A.3J) of the Appendix. 

The normalized electron energy spectrum is shown in Fig. |J2t) ■ In the limit A — ► we have for the 
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Figure 2: Normalized electron energy spectrum (Fr+l) 1 ^ or V° ~ 0-189 (E^o = lOMeV) and a) 

A = m e /m M w 1/207, b) A -> 0. 



x-spectrum 



Jim (Fr+l)- 1 ffi 



= {(1 -yo-x) [-5 - 5 yo + 4y 2 + x (-17 + 14 y ) + 34 x 2 } 
+ 12 x 2 (-3 + 2 x) log(i— ^)} 



(1 - Vo) (7 - 3y + 3y 2 - yo 3 ) + 6 log(y ) 



This is shown in Fig.lJSJ)). 

In the limit y — > we reproduce the spectrum for non-radiative muon decay 



(HP \ non— rad 
— ) =2x 2 (3-2x). 
da; / 



(18) 



(19) 



4 Spectra for Right- and Left-handed Electrons 

The spectra when the final electron is polarized in the state R or L are obtained using Eq.JSJ). We will 
emphasize the case R, which is the one which behaves in an unexpected way in the A — > limit ^ |U] . The 
corresponding spectra for the state L are obtained by subtraction from the unpolarized result L+R discussed 
in Sec|3 
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4.1 Dalitz plot for R-Electrons 

(■J^r) = r " -^T \ ~ {36.4 V +.4(36 y 2 - 60xy 2 - 24 y 3 ) + 12 AX 4 (-8 + y 2 ) 
\dxdy/K 12 it A y 

+ A 2 [ - 288.4 + 96 .4 2 + 192 Ax - 72 A 2 x + 24 Ax 2 + 192 Ay 

- 72 A 2 y + 72 A x y - 18 A 2 x y + 6 A x 2 y + 48 A y 2 + 4 A 2 y 2 

- 52Axy 2 + 9A 2 xy 2 -3 Ax 2 y 2 - 24 Ay 3 - 4 A 2 y 3 + 16 Axy 3 ] 
+ A (A - x) [ - 72 x + 48 x 2 - 72 y + 78 x y + 24 x 2 y + 33 x y 2 

- 10 x 2 y 2 + 12 y 3 - 6 x y 3 + 4 x 2 y 3 ] 
+ C [ - 12 {A - x) (-3 + 2 x + 2 y) (2 x 2 + 2 x y + y 2 ) 
+ A 4 (96.x- 192t/ + 72^y- 24xy + 48y 2 ) 
+ A 2 (288 x - 96 A x - 192 x 2 + 72 A x 2 - 24 x 3 + 72 A y - 264 x y 
- 24 A x y + 72 x 2 y - 48 y 2 - 60 A y 2 + 84 x y 2 + 48 y 3 )] } . 



(20) 



In the limit A — » we obtain 



a^o \dxdy/R 



= r -y(3-2x-2y). 



(21) 



This limiting result coincides with what one obtains on the basis of the helicity-flip splitting function. 
Denoting the electron spectrum in the non-radiative decay (dr/d:ro) non_rad , the spectrum after helicity-flip 
bremsstrahlung is 

( ^ \ / dr \ non_rad ^ 22 

V dxo dz / R V dxo ' h ' 



Noting that y — xq z, x — xq (1 — z), and recalling that 

dr \ non — rad 



we obtain 



exactly as in Eq. l|21|l . 



(—) 

\dxo 1 

(—) =-v 

Vdxdy/R xq Vdxodz/R 



r 2a; 2 (3-2.T ), D M (z) 

dr 



2tt 



r -y(3-2a;-2y), 

7T 



(23) 
(24) 



4.2 Photon-Spectra for R- and L-electrons 

Of particular interest is the photon energy spectrum associated with the right-handed polarization state of 
the electron, which is naively forbidden in the A — > limit. The corresponding spectrum for L-electrons may 
be obtained by subtraction from Ea. (|A.ip . The exact result for (dr/dy)R, is stated in Ea. i|A.5|) . In the limit 
m e/ m /i — > one obtains 

^ r -y(l-y)(2-y), (25) 



a^o \ dy ) R 



which is the result derived in [l] on the basis of the helicity-flip fragmentation function D^z). The spectrum 
for A 7^ is plotted in Fig.©. The striking feature is the appearance of a divergence at y = 0. This is 
connected with the fact that for m e ^ 0, helicity is not a good quantum number. As a consequence, 
ordinary (non-flip) bremsstrahlung, diverging as 1/y at small y, is partly transmitted to the final state with 
electron polarization s e — sf. This contribution vanishes in the limit A — > 0, leaving behind the anomalous 
(helicity-flip) contribution in Eq. (|25|l . This anomalous bremsstrahlung integrates to a decay width 



r R = r — (1 - yo) (1 + 2/o - 3yo 2 + yo A ) , 

4 7T 



(26) 



which is the contribution of right-handed electrons to the muon decay width in the A — * limit 1 J . The exact 
expression for Tr^/o) and the corresponding left-handed rate Tl(?/o) are plotted in Fig. 01. From these, we 
obtain the longitudinal polarization 
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Figure 3: Photon energy spectrum (^p) R for right-handed electrons in units of Tq for a) A = 1/207, b) 
A • I). 
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Figure 4: Integrated rates for R- and L-electrons Fr i l(?/o) m units of Tq as function of minimum photon 
energy y . 



P L (yo) 



Tr(m) - r L (y ) 

^r(vq) + r L (2/o) 



(27) 



plotted in Fig.JSJ). This figure demonstrates the presence of a significant right-handed electron component 
in radiative muon decay, which becomes as probable as the left-handed component at high photon energies. 
To illustrate the approach to the chiral limit A — > 0, we have also shown in Fig.© the polarization of the 
electron in the decay t~ — > e~u e p T j. In the chiral limit A — > the longitudinal polarization goes from 
Pl = — 1 to Pl = like a step- function at yo = 1. 



4.3 Electron Energy Spectra for R- and L-electrons 

The formula for the electron energy spectrum (dr/da;)^ for right-handed electrons is given in Ea. (|A.7|l . and 
may be compared with the unpolarized spectrum (dr/da;)ri,+L m Ea. (|A.3p We have verified that in the limit 
A — > 0, the right-handed spectrum takes the form 



/dT\ a-^ n a _ ^ _ + 5 ^ _ 2 + x _ 2yo) + 2 ^ . 
V da; / R d 7r 

This is exactly the result derived in Ref.pP on the basis of the fragmentation function Dhf(z)- 



(28) 
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Figure 5: Longitudinal polarization of electron Pj, as function of minimum photon energy yo'. a) [i — > e 
decay, b) r — > e decay. 



5 Discussion of the Collinear Region 

A question of interest is the behaviour of helicity-flip and helicity-non-flip brcmsstrahlung in the collinear 
region, 9 ~ 0, where 9 is the angle between the photon and electron in the final state. To this end, we have 
calculated exactly the distribution (dr/dcos(0))R,L- The analytic results are, however, too lengthy to be 
written down here. The left- and right-handed spectra for A = 1/207, in the small angle region, are plotted 
in Fig.©. The longitudinal polarization of the electron as a function of 9 is shown in Fig.JTJ, and, changes 




0.9992 0.9994 0.9996 0.9998 1 
cos(0) 

Figure 6: cos(#) spectrum ( d( wg) ) L R Ior A = 1/207 and yo = 0.189 in units of Tq. 



from the V-A value Pl ~ — 1 to the value Pl ~ +1 in the forward direction. From the exact expressions we 
have derived the behaviour in the chiral limit A — > 0. For 9 ^ and A — >• 0, we find 



VdcosfflWi 



6 7T A 



dcos(6*)/R 

dr \ A^o > r a _1 
dcos((9)/L * ° 3^r 2 



r o ^- 77? [-7 + 10 2/o ~ 6 y 2 + 4 y 3 - y 4 - 6 log(y ) 



(29) 
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Figure 7: Longitudinal polarization Pl for A = 1/207 and yo = 0.189. 



Note that, in the chiral limit, the A/6* 3 behaviour of (dr/dcos(0))R refers to the opening angle of the electron 
and photon in the final state. The corrsponding X 2 /6 4 behaviour of the helicity-flip scattering cross section 
(Eq.([Q)) refers to the angle between the photon and the incident electron. 

We have also derived the value of (dr/dcos(#)) in the forward direction (8 — 0), for small A: 



Vdcos(0)/R 

(-^-) 

Vdcos(6>)/L 



9=0 



r ° ? 3^0 h (1 ~ ya)4 (4 + 16 yo ~ 5 y ° 2) 



9=0 



= r ^ T ^(l-yo) 4 (l + 4 2 /o) 



(30) 



We have found that the angular distribution (dr/dcos(0))R connected with helicity-flip bremsstrahlung can 
be simulated by the following ansatz for the differential decay spectrum (where the label "APP" connotes 
"approximate" ) : 



dT 



APP 



2a 



1 



*dx dy dcos(0)/R 



= To Z_x 2 y(3-2x-2y)- 

n [x- cos(0) A) 



(31) 



When integrated over cos(0), this spectrum, in the limit A — > 0, yields the Dalitz pair distribution given in 
Ea. (|21|l . After integration over x and y we obtain an angular distribution (dr/d cos(0))p PP which agrees 
with the exact calculation over a wide range of angles. This formula is particularly simple for yo = 0: 



/ dJtj/o^O) \ A1 
V dcosf6») ) 



dcos(6*) / r 



A^o „ a 5 

— * r - A 

7T D 



1 + (tt 



cos(0) 



sm{6) 



sm\6) 



(32) 



which agrees with the small angle result given in Ea. (|29|) . Similarly, in the forward direction (9 = 0), for 
small values of A, the approximate representation Ea. (|3"T|) reproduces the result given in Ea. ll!IU|) . 



6 Conclusions 

The purpose of this paper was to examine the dependence of the radiative decay /i - — » e^VgV^ on the 
helicity of the final electron, in order to demonstrate the presence of right-handed electrons, even in the 
chiral limit m e — > 0. This non-intuitive feature is a consequence of helicity-flip bremsstrahlung in QED, 
and was discussed in with the help of the helicity-flip splitting function derived in We have now 
calculated the complete helicity-dependent spectrum of this decay from first principles, without neglecting 
the electron mass, and have shown that the results of are reproduced in the limit A — > 0. The presence 
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of right-handed electrons manifests itself in the longitudinal polarization Pi = (Tr — rL,)/(rn + Tl) which 
deviates from the naive V-A value Pi — — 1 even in the limit A — > 0. The dependence of the polarization 
on the photon energy cut yo, Fig.©, shows that right-handed electrons occur predominantly with high 
energy photons. A characteristic signature of wrong-helicity electrons is also present in the opening-angle 
distribution dr/dcos(0), where the collinear region 9 < A/(l — yo) is dominated by right-handed electrons. 
The fact that helicity-flip bremsstrahlung in the chiral limit m e — > can be exactly reproduced by the 
simple and universal splitting function D^{(z) hi many different processes, is suggestive of a relationship to 
the chiral anomaly in QED [HI El El- 



Appendix A 



(-) 

\ dy ) r+l 



2 A 



( da; dy ) 



dr 

dy ) r+l 



r ° Sr -\[(- 102 + 46y - 95y 2 + 109 y 3 - 45y 4 + 21 y 5 - 6y 6 ) (1 - y) 
36 7r y I 

+ 6 A 2 (64 - 250 y + 386 y 2 - 240 y 3 + 67 y 4 - 20 y 5 + 5 y 6 ) (1 - y)- 1 
+ 18 A 4 y (18 + 12 y - 40 y 2 + 11 y 3 + 3 y 4 ) (1 - y)" 2 
+ 2 A 6 (-192 + 622 y - 856 y 2 + 552 y 3 - 105 y 4 - 21 y 5 ) (1 - y)~ 4 
+ 3 A 8 (34 - 64 y - 13 y 2 + 4 y 3 ) (1 - y)~ 4 ] 

+ [24(-3 + 2y-4 2/ 2 + 2y 3 )(l-y) + 72A 2 y(-13 + 192/-7 2 / 2 )(l- 2 /)- 1 



+ 72 A 8 (-1 + 2 y) (1 - y)- 4 ] log (^-) } 



/d£ N 



a 1 



. -j =r -- M - 4 ( u 2 -A 2 ) 

V da; /r+l 7T ^4 



72 A 4 (12 - 27y + 18y 2 - 5 y 3 ) (1 - y)" 1 + 24 A 6 y (5 + 15 y - 18 y 2 ) (1 - y)" 3 

A 

V 

(A.I) 
, (A.2) 



A- 1 + [j)k 2 + log (^"-) h + log Q) log (i— - ) /,-, 



where u — (x + A)/2, and 



36 



(u 2 - A 2 ) (1 - y - u) u 2 ( - 300 + 399 u + 71 u 2 + 2 u 3 + 8 u 4 + 132 y + 63 u y - 10 u 2 y 

- 8 u 3 y + 24 y 2 - 12 u y 2 + 8 u 2 y 2 ) 
+ A 2 it (555 - 259 u - 61 u 2 - 19 u 3 + 87 y - 28 it y + 19 u 2 y 
- 12y 2 + 8uy 2 ) 



A 4 (122 - 67 u - 19 u 2 - 22 y + 19 u y 



fe 3 

fe4 



— it (1 — yo — w) u 2 (5 + 17 u — 34 u 2 + 5 yo — 14 u yo — 4yo 2 ) + 3 A 2 u (6 — 19 it + u 2 — 6 yo — 5 wyo) 

+ 12A 4 (-4-m + 3u 2 ) , 

4 (A 2 -u 2 ) [u 3 (3 — 2 u) + A 2 u (3 - 8 u + 3 u 2 ) + A 4 (-2 + 3u)] , 
4(A 2 +u 2 ) [u 3 (3-2u) + A 2 u(3-8u + 3u 2 ) + A 4 (-2 + 3u)] , 

(A.3) 
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( dy ) r J 2 \ ^ ( dx dy ) r 



-r a 

— 1 o — 

7T 



y (1 - y) (2-y) + h 1 + log (-^) ^ 2 + log 2 (—^) h 3 



where 



hx := 4 A y (-3 + 2 y) + — A 2 (-248 + 176 y + 489 y 2 - 304 y 3 - 53 y 4 - 6 y 5 ) 

18 y 1 - y 

+ | A 3 - (48 - lOy + 15 y 2 + 3y 3 ) + \ A 4 - 1 (-150 + 320y - 209 y 2 + 16y 3 + 29 y 4 - 2y 5 ) 
9 y 4 y (1 — y)^ 

+ 1 A 5 - (32 - 12 y - 3 y 2 ) + 4 A 6 - , n 1 .„ (-120 + 412 y - 465 y 2 + 114y 3 + 18 y 4 ) 
9 y 18 y (1 - y) 3 

+ 4 A " " 7T^4 ( 38 - 140 y - 21 y 2 + 6 y 3 ) , 
36 y (1 - y) 4 

^2 := | A 2 - -J— (8 + 13y - 30y 2 + 7y 3 + 5y 4 ) + A 4 - — !— (-16 + 22y - 19y 2 + 3y 3 ) 
3 y 1 - y y 1 - y 

+^ 6 7T^(- l - 6 ^^ 2 )+^ 8 ^7J^ F (l-4 ^ /), 

ft 3 := 2 A 2 (1 - y) (-3 - y) + 2 A 4 ± (-2 - y + y 2 ) . 



V dx / R 7T 



( ,: - j x + log j 2 + log (- — -) j 3 + log log (- 1 



yo / V A/ v y 



J4 



(A.5) 
(A.6) 



where 

1 1 

3d u 



3 (1 - y - u) {u 3 [30 + 12 u 2 + u (-42 - 12 y ) + 30 y - 24 y 2 ] 

+ A 2 m 2 [264 - 168 y - 24 y 2 + it (261 - 15 y - 12 y 2 ) 

+ u 2 (-155 - 26 y + 4 y 2 ) + u 3 (-5 + 5 y ) + u A (-5)] 
+ A 4 u [ - 555 - 87 y + 12 y 2 + u (199 - 32 y + 4 y 2 ) 

+ w 2 (73 + 5y ) + u 3 (-5)] 
+ A 6 [- 122 + 22yo-8yo 2 + 10u(4-y ) + 10 U 2 ]}, 

2 r 

7a2 {w 2 (1 - yo - m) [ - 5 - 5 y + 4 y 2 + u (-19 + 12 y + 4y 2 ) + w 2 (-13 - y ) + 13 1 
+ 6 A 2 u (1 - y - u) [ - 3 + 3 y + u (2 + 3 y ) + u 2 (-4 + y ) + it 3 ] 
+ 24 A 4 (2-u)(l + u)(l-ito-u)}, 



j 3 := 4 ^ A 2 (1 - u) (A 2 - u 2 ) (-2 A 2 + 3 u - u 2 ) , 



j 4 := 4-^A 2 (l- U )(A 2 + u 2 )(-2A 2 + 3u-w 2 ). 



(A.7) 
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